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ABSTRACT 

TeV-blazars are known as prominent non-thermal emitters across the entire 
electromagnetic spectrum with their photon power peaking in the X-ray and TeV- 
band. If distant, absorption of 7-ray photons by the extragalactic background 
light (EBL) alters the intrinsic TeV spectral shape, thereby affecting the over- 
all interpretation. Suzaku observations for two of the more distant TeV-blazars 
known to date, 1ES 1101-232 and 1ES 1553+113, were carried out in May and 
July 2006, respectively, including a quasi- simultaneous coverage with the state 
of the art Cherenkov telescope facilities. We report on the resulting data sets 
with emphasis on the X-ray band, and set into context to their historical be- 
haviour. During our campaign, we did not detect any significant X-ray or 7-ray 
variability. 1ES 1101-232 was found in a quiescent state with the lowest X-ray 
flux ever measured. The combined XIS and HXD PIN data for 1ES 1101-232 
and 1ES 1553+113 clearly indicate spectral curvature up to the highest hard 
X-ray data point (~ 30 keV), manifesting as softening with increasing energy. 
We describe this spectral shape by either a broken power law or a log-parabolic 
fit with equal statistical goodness of fits. The combined 1ES 1553+113 very high 



1 W.W. Hansen Experimental Physics Laboratory & Kavli Institute for Particle Astrophysics & 
Cosmology, Stanford University, 452 Lomita Mall, Stanford, CA 94305, USA; afr@stanford.edu, 
Luigi.Costamante@stanford.edu, olr@stanford.edu 

2 Kavli Institute for Particle Astrophysics and Cosmology, Stanford University, Stanford, CA 94305, 
and Stanford Linear Accelerator Center, 2575 Sand Hill Road, Menlo Park, CA 94025, USA; made- 
jski@slac.stanford.edu 

3 Universitat Wurzburg, D-97074 Wurzburg, Germany; dorner@astro.uni-wuerzburg.de 



-2- 



energy spectrum (90 — 500 GeV) did not show any significant changes with re- 
spect to earlier observations. The resulting contemporaneous broadband spectral 
energy distributions of both TeV-blazars are discussed in view of implications for 
intrinsic blazar parameter values, taking into account the 7-ray absorption in the 
EBL. 

Subject headings: X-rays: galaxies - galaxies: active - BL Lacertae objects: indi- 
vidual (1ES 1101-232, 1ES 1553+113) 



Introduction 



Blazars are among the most extreme sources in the high energy sky. They constitute a 
subclass of the jet-dominated radio- loud population of active galactic nuclei (AGN), which 
in turn are sub-divided into flat spectrum radio quasars (FSRQ) and BL Lac objects where 
the absence or depression of strong emission lines characterizes the latter. Their observa- 
tional properties include non-thermal continuum emission and irregular variability across 
the whole electromagnetic band, ranging from subhours at gamma-ray energies to week-long 
time scales in the radio band, often high optical and radio polarization, and a core-dominated 
radio morphology. In some sources, superluminal motion has b een detected, indicating a rel- 
ativistically enhanced emission region close to the line-of-sight ( jBlandford fc Rees 19781 ). As 
a consequence of this, any jet emission is highly beamed. 

The spectral energy distribution (SED) of blazars shows two broad peaks in the vF v 
representation. The low energy hump is generally agreed to stem from synchrotron ra- 
diation of relativistic electrons and/or positrons in the jet whereas the origin of the high 
energy one is still under debate. Leptonic models explain the complete non-thermal SED 
as synchrotron and inverse Compton emission from relativistic electrons or pairs which 
upscatter their self-produced synchrotron radiation or external photon fields. Hadronic 
models produce the high energy peak either via interac tions of relativis tic protons with 
matter (Pohl fc Schlickeiser 2000 ). ambient photons (e.g. Mannheim 19931), magnetic fields 



jAharonian 200oh and/or both, magnetic fields and photons (e.g. iMiicke &: Protheroe 2001 



Miicke et al. 20031 ) 



Almost all of the 19 blazars (except for BL Lacertae ([Albert et al. 2007bl ) and 3C 279 
(ITeshima et al. 20071 )) detected reliably to date at the TeV energy belong to the subclass of 
high-frequency peaked BL Lacs (HBLs), with their two vF v peaks typically at UV/X-rays 
and TeV-energies. Until recently blazars have been detected preferentially in their flare state, 
owing to instrumental limitations. With the contemporary Cherenkov telescopes (such as 
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H.E.S.S., MAGIC, VERITAS) it is possible to study also low activity states (subsequently 
referred to "quiet state"). 

TeV-blazars have also been used successfully as a powerful tool to probe the extragalactic 
background light (EBL) at IR/optical energies via the absorption of 7-rays in the cosmic 
diffuse radiation field. The most stringent constraints on the EBL density at ~ 2/im stem 
from observations with the H.E.S.S. telescope syst em of the high-redshift, hard-spectrum 



TeV-blazar 1ES 1101-232 (lAharonian et al. 2006al ). The most distant (redshift z > 0.1 



TeV-blazars established to date are 1ES 1101-232 (z = 0.186), 1ES 1218+304 (z = 0.182), 
1ES 1011+496 (z = 0.212), 1ES 0347-121 (z = 0.185), 1ES 0229+200 (z = 0.140), H 2356- 
309 = 0.165), H 1426+428 (z = 0.129), PKS 2155-304 (z = 0.116), 1E S 1553+113 



z > 0.25: IPerlman et al. 19961 : iFalomo fc Treves 19901 : lAharonian et al. 2006bf ) 



Extreme HBLs, as is the case for most very high energy (VHE) blazars, are character- 
ized by their spectral extension to extremely high particle energies. TeV-blazars have been 
regularly observed in the past by X-ray instruments, however, not many have been detected 
so far in the hard X-ray range where the highest energy particles leave their imprints. The 
sensitivity reached by the hard X-ray detector HXD onboard the Suzaku satellite allowed for 
the first time to study in detail possible new spectral features and curvature at > 10 keV, 
even in a blazar quiet state, which may help understanding the properties and origin of the 
VHE emission. 

We report here on our ~ 50 ksec Suzaku observations on each of two more distant TeV- 
blazars known to date: 1ES 1101-232 and 1ES 1553+113. These observations were covered 
by quasi-simultaneous observations with the H.E.S.S. telescope system and the MAGIC 
telesc ope. None of these sources have shown in the pa st significant variability in the TeV- 



band (jAharonian et al. 2006al : lAharonian et al. 2006bl ). likely owing to their low flux close 



to the detection limit of the TeV facilities. 

The paper is organized as follows: We provide a brief overview of both TeV-blazars in 
Sect. 2. Sect. 3 is devoted to a technical description of our Suzaku data analysis and its 
results. In Sect. 4 we report on the quasi-simultanuous multiwavelength coverage during the 
Suzaku observations of 1ES 1101-232 and 1ES 1553+113. In Sect. 5 we discuss our results 
in particular in context to their historical behaviour. Sect. 6 summarizes the main results of 
this work. 
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2. The VHE blazars 1ES 1101-232 and 1ES 1553+113 



The X-ray selected BL Lac object (XBL) 1ES 1101-232, hosted by an elliptical galaxy at 
redshift z = 0.186 with estimated magnitude uir — 16.41, one of the brightest BL Lac hosts 
known so far, is presumably located in a galaxy cluster. It has been detected by various X-ray 
instruments: the Ariel-5 X-ray satellite discovered this source, HEAO-1, Einstein, ROSAT, 
BeppoSAX, RXTE, XMM, SWIFT, Suzaku (this paper) j ust to name a few. Its st rong X-ray- 
to-optical flux ratio led to the classification of an XBL dScaroa fc Falomo 1997h . which has 
later been revised t o a high-frequency p eaked BL Lac object (IDonato et al. 20011 ) within the 

HBL-LBL scheme (iGiommi et al. 19951 ). The radio maps of 1ES 1 101-232 show a one-sided, 

not well-collimated jet structure at a few kpc distance from the core (ILaurent-Muehleisen et al. 19931) 



The optical core emission is typically varying on time scal es of months (e.g. 



Remillard et al. 198 



J, 



although occasional intraday flares have been reported (IRomero et al. 19991 ) . Based on the 
typical broadband char acteristics of VHE AGN several authors have predi cted this source to 



be a TeV-blazar (e.g. IWolter et al. 20001 ; ICostamante &: Ghisellini 20021). Indeed, it wa s 



recently detected by the H.E.S.S. Cherenkov telescope system ( lAharonian et al. 2006bl ). 
At MeV-GeV en ergies, however, EGRET reported only an upper limit from 1ES 1101-232 
Lin et al. 1996h . 



Owing to the non-detection of any lines nor resolving its host galaxy the redshift of 
the XBL/HBL 1ES 1553+113 is to date still un known. Redshift limits , though, have been 
imposed either through HST results (z > 0.78: ISbarufatti et al. 20051 ) on the assumption 
that the nucleus' host is a typical BL Lac host, or t hrough arguments involy ing TeV pho- 
ton absorption in th e EBL ((z > 0.25, z < 0.74: lAharonian et al. 2006bl ): (z < 0.42: 
Mazin fc Gobel 20071 )) and no anomalies in the source intrinsic high energy SED. Like 1ES 
1101-232, 1ES 1553+113 is a strong X-ray source and has been visited by several X-ray 
instruments, however, focusing mainly towards the soft part of its X-ray spectrum (Einstein, 
ROSAT, BeppoSAX, R XTE, XMM). Its weak radio flux (a few tens of a mJy at 5GHz: 
Falomo fc Treves 199ph among the blazar class l eads to its classification of an XBL and 
HBL ( jFalomo fc Treves 1990t IDonato et al. 20051 ). Results from a 3-weeks multifrequency 
campaign in 2003 covering the radio-to- X-ray range (lOsterman et al. 20061 ) where a strong 
X-ray flare was detected, indicate an extremely high synchrotron cutoff energy, superceding 
in this respect even well-known extreme HBLs like Mkn 421, Mkn 501 or PKS 2155-304. Thi s 
places 1ES 1553+113 right at the frontline of the class of extreme HBLs ( IRector et al. 20031 ). 
Despite its unknown redshift, it was considered as a promising source for detection at TeV 
energies based on its spectral broad band properties. The H.E.S.S. collaboration reported 
recently it s discovery at VHEs with a integral flux detection on a 3.8cr level, and a very sof t 
spectrum ( lAharonian et al. 2006bt). This was con firmed by MAGIC ( lAlbert et al. 2007al ). 
The Third EGRET catalog (IHartman et al. 19991 ) does not list this blazar as identified in 
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the MeV-GeV energy band. 



3. Suzaku observations and data analysis 

The two object s have been observed by Suzaku's X-ray Imaging Spectrom eter (XIS: 



Koyama et al. 20071 ) and the Hard X-ray Detector (HXD: iTakahashi et al. 20071 ) in May 25- 
27 2006 for a total of 53.1 sec on-source for 1ES 1101-232, and in July 24-25 2006 for a total of 
41.1 sec on-source for 1ES 1553+113. XIS is a CCD instrument with three front-illuminated 
and one back-illuminat ed CCD camera at the focal plane of the four X-ray telescopes (XRT: 



Serlemitsos et al. 20071 ). and sensitive in the energy range 0.2-12 keV. During our observa- 
tions of the two TeV-blazars, XIS was operated in the 3x3 and 5x5 observation modes. 
No charge injection procedures were carried out during these observations. 

The HXD consists of an array of 4 x 4 detectors (well units) and 20 surrounding crystal 
scintillators for active shielding. Each unit consists of four Si PIN diodes (PIN), sensitive in 
the 10-70 keV energy range, and four GSO/BGO phoswitch counters (GSO) for detecting 
photons in the 40 - 600 keV regime. Since the focus of these observations were the hard part 
of the X-ray energy regime, the HXD was chosen as the nominal point to ensure optimal 
sensitivity. 

One day before our Suzaku observations of 1ES 1101-232 started, the HXD team started 
changing the bias voltage, affecting potentially HXD's UNIDIDs to 3, to suppress the rapid 
increase of noise events noted at this time. In particular, on May 26 2006, the bias voltage 
of the HV-P0 was increased in two steps, affecting the measurements of the PIN during the 
1ES 1101-232 observations: we found two huge background 'flares' in UNIDID 1, towards 
the end of this observation beyond MJD 53941.7. 

Our data analysis is based on the version 2 processed data with the standard mkf 
filtering, and using tools of the HEASoft version 6.4. The CALDB release from February 1 
2008 contains the calibration files used in this analysis. For the data reduction we followed 
the "The Suzaku Data Reduction Guide"0 and specifically Ver. 2.0 if not noted otherwise. 



1 http:/ /suzaku. gsfc.nasa.gov/docs/suzaku/analysis/abc/ 
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3.1. XIS data reduction and analysis 

We started the XIS data analysis with the unscreened version 2 of the data. We com- 
puted the XIS pulse invariant and grade values and updated the CTI calibration via the tool 
xispi. The subsequent reduction, including cleaning of hot and flickering pixels, selection of 
grades 0, 2, 3, 4, and 6, filtering out the epochs of high background (e.g. due to the satellite 
crossing through the South Atlantic Anomaly) was performed using the tool xselect (and 
specifically a script xisrepro. Source events were then extracted in a circle around the 
source with radius of ~ 260". The extraction of background events were performed from a 
region devoid of any obvious X-ray sources. We ensured that the use of other background re- 
gions within the XIS field does not significantly affect the results from our spectral analysis. 
We generated the corresponding response (rmf ) and auxiliary (arf) files by utilizing the tools 
xisresp and xissimarf gen, using the default energy steps; the latter is a Monte-Carlo code 
provided by the Suzaku team, and it also accounts for the hydrocarbon contamination that 
has been built-up in the optical path of each XIS sensor. With the most recent calibration 
the spectra do not show calibration anomalies down to ~ 0.4 keV, which we therefore in- 
cluded into the analysis. Our spectral results do not change significantly even when limiting 
the XIS bandpass to above 0.7 keV. 



3.2. HXD data reduction and analysis 

As was the case for the XIS, here we also used the version 2 processed data. Here, we 
started with the filtered data, where standard cleaning criteria were used; however, for the 
extraction of the spectrum of 1ES 1101-232 we selected only the data before MJD 53941.7 to 
avoid any contamination with the two huge background flares. We also extracted the non- 
X-ray background spectrum from the background files provided by HEASARC, but applied 
the same time selection criteria as to the source files. The source spectrum was dead-time 
corrected using the tool hxddtcorr; we note that the application of hxddtcorr resulted in a 
live time of 93.4% for the 1ES 1101-232 and 93.8% for the 1ES 1553+113 observations. We 
applied analogous procedure to the GSO data, but since neither source was detected in the 
GSO, we did not use those data in the subsequent analysis. 



For the spectral analysis of the combined t ime- integrated XIS and PIN data we have 
taken into account the cosmic X-ray background (IGruber et al. 19991) n ormalized to the PIN 
field-of-view, and a XIS/PIN normalization of 1.12 (jlshida et al. 20071 ). 
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3.3. Results 

We detected 1ES 1101-232 and 1ES 1553+113 with Suzaku up to energies above ~30 keV. 
This is the highest X-ray energy measurement so far obtained for 1ES 1553+113. 

The combined XIS-PIN data has been fitted to both, simple and broken power laws. For 
both sources we found a significantly improved x 2 -vahie when fitting to a broken power law as 
compared to a single power law. The absorbing column densities were fixed to the respective 
Galactic values here (N H = 5.76 • 10 20 cm" 2 for 1ES 1101-232, N H = 3.67 • 10 20 cm" 2 for 
1ES 1553+113). The best-fit parameters and goodness-of-fits are summarized in Table [TJ 
Both sources show the break to be around 1.4 keV with 1ES 1553+113 having a somewhat 
steeper overall spectrum than 1ES 1101-232. We determined the average 2-10 keV flux from 
1ES 1101-232 to be 1.68 ■ 10~ n erg cni'V 1 , and the 10-30 keV flux to be 7.4 ■ 10~ 12 erg 
cm _2 s _1 (source flux only). This is the lowest flux level so far measured from this source. 
The average 2-10 keV source flux from 1ES 1553+113 is 3.54 ■ 10~ n erg cm~ 2 s _1 , about a 
factor 2 higher than for 1ES 1101-232. The 10-30 keV flux is at 1.35 ■ 10~ n erg cm" 2 s^ 1 . 
This X-ray flux level is comparable to the 2001 XMM data, and lies at an intermediate value 
between the so far highest (found during the 2005 SWIFT observations) and lowest (detected 
during the 2003 RXTE campaign) flux state of 1ES 1553+113. The best-fit and residuals 
for both sources are shown in Fig. [TJ and [3J 

We tested for variability in both Suzaku data sets using the y 2 method. No flux vari- 
ability on subhour or intraday time scales were found (Fig. [2], |4|). The lack of variability, 
otherwise often observed on these time scales during high flux states in TeV-blazars such 
as Mkn 501, Mkn 421, PKS 2155-304, etc., may suggest the sources to have lingered at a 
quiescent activity level in 2006. 

Log-parabolic shapes have bee n proposed to provide a better description of X-ray spec- 



tral shapes in HBL-type AGN (e.g. iPerlman et al. 20051 ). A x 2 -test to the shape dN/dE oc 



_g;[-r-/3iog(£)] Q f com bined XIS-PIN spectrum revealed a statistically comparable good- 
ness of fit as compared to the broken power law representation up to the highest reported 
energy point for both sources. Best-fit parameters for the log-parabolic shapes are in Ta- 
ble [TJ Although the best-fit curvature value for both sources is similar, we note a ten- 
dency of a further steepening in the PIN data with respect to the model (log-parabolic 
or broken power law shape) describing the 1ES 1553+113 data set (see e.g. Fig. [3]). A 
log-parabolic shape spanning the whole dec line of the synchrotron component can be the re - 



sult of e nergy- dependent acce leration (e.g. iProtheroe fc Stanev 19991; iMassaro et al. 2 004) 



episodic (IPerlman et al. 20051 ) or fluctuating particle energy gains (ITramacere et al. 20071 ). 
If this interpretation holds, the present result indicates acceleration operating also in a quiet 
state of source activity. According to Protheroe & Stanev (1999) the spectral shape of 
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1ES 1101-232 is in agreement with a cutoff due to either E 2 particle energy losses or a finite 
size of the emission region for an energy dependence of acceleration and/or escape oc E^ 1 ^ 
or weaker. 

Our investigation of the Suzaku data found no convincing signatures of spectral harden- 
ing up to the highest detected X-ray energy point, as would be expected from a synchrotron 
emitting Compton loss dominated source in an energy range affected by the Klein-Nishina 
decline of the cross section, or the onset of the 7-ray component. We therefore conclude that 
the measured X-ray radiation is entirely due to the synchrotron component, and that losses 
are either not Compton loss dominated in the X-ray band, or are still in the Thomson loss 
regime. 
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Fig. 1. — Data and folded broken power law model of 1ES 1101-232 that fits the data best 
( x 2 red = l.H) and residuals. Parameters: I\ = 2.04 ± 0.02, T 2 = 2.32 ± 0.02, E hreak = 
1.37 ± 0.08keV. The four curves below 10 keV correspond to the four XIS detectors while 
the curve above 10 keV denote the PIN data. 
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Fig. 2. — Background corrected lightcurve of XIS data sampled with 256 sec binning. A \ 2 
fit to a straight line indicates no signatures of variability of the source during the Suzaku 
observations. 
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Fig. 3. — Data and folded best fit (x? e d = 1-08) broken power law model of 1ES 1553+113. 
Parameters: T x = 2.13 ±0.02, T 2 = 2.42 ±0.01, £ break = 1.38 ±0.06 keV. Note that the PIN 
data lie systematically below the model indicating a further steepening. For a description of 
the individual curves see Fig. [TJ 
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channel energy (keV) 

Fig. 4. — Background corrected XIS lightcurve sampled with 256 sec binning. A %2-test 
indicates a slight preference of representing the XIS lightcurve with a linear fit better than 
with a constant flux. We found no signatures of short time variability during the Suzaku 
observations. 
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Object 


Ti 


r 2 


-E'break 


-^red 


r 


a 


-^red 








keV 


(d.o.f.) 






(d.o.f.) 


1101-232 


2.04 ±0.02 


2.32 ±0.02 


1.37 ±0.08 


1.11 (3328) 


2.10 ±0.01 


0.26 ±0.02 


1.10 (3329) 


1553+113 


2.13 ±0.02 


2.42 ±0.01 


1.38 ±0.06 


1.08 (4354) 


2.19 ±0.01 


0.26 ±0.01 


1.06 (4355) 



Table 1: Spectral fit results from the combined XIS-PIN Suzaku data of a broken power law 
fit with photon indices Y\, r 2 and break energy -Ebreak and of a log-parabolic shape with 
photon index T and curvature (5. Galactic column densities has been used for these fits. 
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4. Simultaneous multifrequency observations 
4.1. Observations in the optical band 

The KVA observatory on La Palma is monitoring 1ES 1553+113 reg ularly since early 



2006. KVA observed this source between 21 and 27 July 2006 in the R-band (I Albert et al. 20081 ) 
The KVA pointing date that falls into the Suzaku observing time is MJD=53941. 43495. 
The measur ed flux was correcte d for Galactic extinction using a R-band extinction of A = 



0.172 mag (jSchlegel et al. 19981 ). Since a host galaxy contribution for 1ES 1553+113 is 
negligible, no host galaxy correction for the flux has been made. A flux value of S op t = 
14.28 ± 0.23 mJy was derived. 



4.2. Simultaneous TeV-obser vat ions 



1ES 1101-232 was visited by the High Energy Stereoscopic Syst em (H.E.S.S.) atmo- 
spheric Cherenkov telescope array several tim es between 2004 and 2007 (jAharonian et al. 2007a 
Aharonian et al. 2006al ; iBenbow et al. 20071 ). In May 2006 H.E.S.S. observations were sched- 
uled to cover the epoch of the Suzaku measurements towards 1ES 1101-232. The 4.3 h 
good-quality data revealed only a marginal detection of 2.9cr (51 e xcess events) and a 
corresponding flux of J (> 260GeV) = 3.2 ± 1.4 stat x lO^cm-V 1 faenbow et al. 20071 : 
Aharonian et al. 2008bf ) . The complete 2006 data set showed a 3.6<r excess in 13.7 h live 
time. In comparison, the 2004-2005 TeV-data on this source revealed an integral flux of 
/(> 200GeV) = 4.5 ± 1.2 x lO^cm-V 1 (10.1a) JAharonian et al. 2007ah . The low X-ray 
flux level during the Suzaku compaign was therefore accompanied with a comparable VHE 
flux when compared to the 2004/05 activity state. All 2004/05 VHE data are compatible 
with a diffe rential power law spectru m dN/dE = (5.63 + 0.89) • 10 _13 (E/TeV) _ ^ 2 ' 94±0 ' 20 ^cm -2 
s" 1 TeV" 1 dAharonian et al. 2007ah . 



The Suzaku observations of 1ES 1553+113 in July 2006 were also accompanied by si- 
multaneous VHE observations, with H.E.S.S. and MAGIC. The total available live time data 
in July 2006 from H.E.S.S. encompasses 10.1 hours, out of which 3.1 hours live time obser - 
vations were carried out simultaneously with the Suzaku pointing (jAharonian et al. 2008al ). 
An excess of 101 events were found, corresponding to a significance of 3.9a. For the total 
July 2006 data the significance incre ased to 6.7a with an inte grated flux of /(> 300GeV) = 



(4.33 ± 0.94) x 10 erg cm 2 s (jAharonian et al. 2008af ). A daily light curve did not 



show any convincing evidence of variability, as did a run-by-run light curve from the Suzaku 
epoch observations (see Fig. [5]). The low significance found by the H.E.S.S. observations 
during the Suzaku campaign did not justify to produce a reliable fit to the 7-ray data 
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( lAharonian et al. 2008al ). The whole July 2006 data set was reported with a power law 
fit with spectral photon index T 7 = 5.0 ± 0.7. The spectral data points derived from this 
epoch are shown in Fig. [6J This is in very good agreement (within 2a) with the TeV-spectra 
derived from the past observations in 2005 a nd 2006. The combined 2005/2006 data set is 
compatible with a spectral index of T 7 ~ 4.9 (lAharonian et al. 2008al ). 



In July 2006 MAGIC was also observing 1ES 1553+113 for a total of 9.5 hours. These 
observations suffered from a sand-dust coming from the Sahara (Calima). The effects of 
this weather condition on the data has been corrected using extinction measurements from 
the Carlsberg Meridian Telescope ([Albert et al. 20081. for det ails). With a live time of 8.5 
hours the signal significance reached 5.0<j ( Albert et al. 20081 ). enabling the compilation of 
a VHE spectrum for this time span. The spectrum above ~ 90 GeV can be fitted by 



dN/dE= (1.4±0.3)10- 10 (£/200GeV) 



-4.1±0.3r 



TeV s cm , and is in good agreement with 



the corresponding H.E.S.S. data points of July 2006. For the 2 days of contemporaneous 
X-ray/VHE coverage the MAGIC data allowed only a ~ 2a result in 2.4 hours live time. 
The spectral data points from these 2 days observations lie within the uncertainties of the 
overall July 2006 spectrum, which therefore will be used in the following as an adequate 
description of the VHE spectrum for the time of the Suzaku observation. 

Fig. [7] shows the spectral data points of the July 2006 MAGIC measurements together 
with the H.E.S.S. data points for this time span, covering the energy range of ~ 90—500 GeV. 
A fit to the combined H.E.S.S. and MAGIC data points during the July 2006 observations 
revealed a simple power-law representation to be adequate. When compared to broken power 
law with exponential cutoff spectral shapes in the 90-500 GeV energy range, the goodness 
of fit does not improve noticeable. The day-by-day light curve does not show any signs of 
variability (see Fig. EJ), although the average flux > 150 GeV lies systematically above the 
corresponding H.E.S.S. > 300 GeV flux (because most of the VHE signal stems from below 
200 — 300 GeV). Also MAGIC confirms the sur prising stability of t he spectral shape of this 
source, even when compared to year 2005 data (I Albert et al. 20081 ). 



Up to now 1ES 1553+113 is lacking any confirmed redshift measure ment. The initially 
proposed redshift of z — 0.36 turned out an mis-interpretation of line data (IFalomo fc Treves 19901 ) 
Const raints have been set only from indir ect methods, e.g. from limits on host galaxy im- 
ages (jUrry et al. 200(1 iTreves et al. 20071 ) one arrives at z > 0.25. Recent spectroscopic 
obser yations improved here to z > 0.09 by deriving limits of the host galaxy absorption 
lines ( jSbarufatti et al. 20061 ). Gamma ray absorption in the EBL followed by pair pro- 
duction in the VHE range offers a further possibility to infer on redshift limits: The re- 
construction of the emitted 7-ray spectrum from the observed one with imposing limits on 
the hardness of the intrinsic source spectrum leads to the sought after redshift constraint. 



Such and similar concepts have been appli ed to past observations, and yielded z < 0.7 4 
flAharonian et al. 2006bl : I Albert et al. 2008^ . z < 0.69 and z < 0.42 Jviazin fc Gobel 2007h . 



In Fig. [7] we show the de-absorbed TeV-spectrum of 1ES 1553+113 for a source redshift 
of z = 0.3 and using various models for the EBL: the P0.45 flux level from Aharonian et al. 
(2006a), the "best-fit" model of Kneiske et al. (2004) and the latest version of the baseline 
model of Stecker et al. (2006). The P0.45 flux level corresponds to a EBL flux density 
with a phenomenological spectral shape that describes the EBL data points satisfactorily on 
the level of the integrated galaxy counts. All models support about the same 2/zm EBL-flux 
density, however, differ in their extension to lower and higher energies. Note that the baseline 
model flux at < 2/im is not based on observational constraints, but rather represents the 
extension of the Stecker et al. (2006) diffuse IR background calculation into the optical/UV 
band, and should therefore be considered with caution in this energy range. The resulting 
power spectra appear all at around photon index ~ 3 when de-absorbed by either model with 
the baseline model providing the highest intrinsic VHE flux level up to ~ 450 GeV. Still all 
considered EBL models lead to comparable photon powers delivered through the high and 
low energy SED humps. Stepping up in redshift the de-absorbed spectra are getting harder 
quickly. At a source redshift of z — 0.67 a spectral photon power law oc E^ 1 is reached 
at > 300 GeV (see Fig. [7]), the asymptotically hardest spectrum from any conservative 
particle acceleration scenario. We therefore consider a redshift of z > 0.7 as unlikely for 1ES 
1553+113 if the EBL flux density is at the lowest possible value (integrated galaxy counts). 
A search for the redshift that corresponds to an intrinsic E~ 2 photon spectrum concludes at 
z = 0.5 (see Fig. UJ). This represents the borderline for the VHE bump to lie at either tens of 
GeVs or >TeV. When considering redshifts < 0.3, the location of the VHE SED peak moves 
to < 90 GeV with an estimated power output at 7-rays again comparable to the synchrotron 
power peak. The intrinsic photon spectral index is then ~ 3.1. 
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Fig. 5. — Combined VHE light curve from H.E.S.S. (tri angles; I(> 300GeV); daily; 
average flux is 5.8 ± 1.7 sta t ± 1.2 sys t x 10 -12 cm -2 s _1 ; ( lAharonian et al. 2008al )) and 
MAGIC (diamonds ; /(> 150GeV); daily; average flux is 27.0 ± 1.7 stat x 10" 12 cm" 2 s _1 ; 
( lAlbert et al. 20081 )) of 1ES 1553+113. No indication of daily variability is apparent. Dashed 
vertical lines correspond to the start and end times of the Suzaku observations. Inlay: 
VHE light curve from H.E.S.S. during 24 and 25 July 2006 (run-by-run; average flux is 
4.6 ±0.6 sta t ±0.9 syst x 10~ 12 cm -2 s _1 ; (lAharonian et al. 2008a] )). No indication of variabil- 
ity is apparent on this time scale. 
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Fig. 6. — Observed VHE spectrum from 2005 (open symbols) H.E.S.S. (circles) and MAGIC 
(squa res) observations of 1ES 1553+113, compar ed to the July 2006 VHE data (filled sym- 
bols; lAharonian et al. 2008al : I Albert et al. 20081 ). 
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Fig. 7. — Combined VHE spectrum from H.E.S.S. (blac k triangles; lAharonian et al. 2008al ) 
and MAGIC (black diamonds; I Albert et al. 2008^ of 1ES 1553+113. The com- 
bined fit (dashed line) follows a power law dN/dE = (2.192 ± 4.0742) x 
10- 23 ( J B Te v/267.912)- 4 - 068±0178 cm- 2 s- 1 TeV- 1 (x 2 ed = 1.50) where E TcY is the energy in 



TeV. Green /light grey points belong to the de- absorbed VHE data usi ng the minimal EBL 
(P0.45 from Aharonian et al. 2006a ). red/dark grey points are for the ( Kneiske et al. 2004h 
"best-fit" EBL-model and blue/medium grey points correspond to the (jStecker et al. 20061 1 
baseline EBL-model used for the de-absorption and assuming a source redshift of z — 0.3. 
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Fig. 8. — The intrinsic (de-absorbed) TeV-spectrum of 1ES 1553+113 for a source redshift 
of z — 0.67 (blue/black), z = 0.5 (green/light grey ) and z = 0.3 (red/dark grey) using the 
minimal EBL (P0.45 from lAharonian et al. 2006al ). The source redshifts at z > 0.3 have 



been chosen to lead to intrinsic E 1 (z = 0.67) and E 



0.5) photon power law spectra. 
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5. Discussion 

5.1. Historical X-ray behaviour of 1ES 1101-232 and 1ES 1553+113 

1ES 1101-232 has been detected during the past 15 years by a range of X-ray instru- 
ments, both covering the soft (ROSAT, XMM) and extending into the hard X-ray range 
(BeppoSAX, RXTE, SWIFT) at various sensitivity levels. Except for the continuous mon - 



itoring by RXTE's ASM and the 19.6 ksec monitoring by XMM (jAharonian et al. 2007al ). 
most observations were carried out in snapshot mode, with each observation encompassing 
typically 1-5 ksec. The continuous Suzaku observation presented here constitutes therefore 
the first long (> 30 ksec) high- sensitivity X-ray data set of this source. None of the past 
observations showed significant short time variability (<hour). Although the continuous 
mode observations of Suzaku generally allow to probe intraday variability, we find no sig- 
nificant variability of this time scale during our observations. Year-by-year variations as 
derived from the historical X-ray light curve are less than a factor ~ 2. The highest ever 
measured X-ray flux ley el has been reported for 2005 from the RXTE-PCA measurements 



(jAharonian et al. 2007al ). There, in eleven consecutive nights the source was observed each 
night for ~ 10 ksec. No hints for strong flux variability were found, and the AGN was catego- 
rized to be in a non-flaring state. If this is the case, 1ES 1101-232 was likely never observed 
in a flaring state during the past 15 years. Only 4 months later SWIFT saw the sou rce in 



a smooth decline of flux level with no sub-hour variations (e.g. iMassaro et al. 20081 ). Our 
Suzaku observations, made approximately 10 months later, found the lowest ever measured 
X-ray flux from 1ES 1101-232. 

So far, essentially all sensitive X-ray observations of 1ES 11 01-232 prefer a curve d or 



broken power-law representation instead a simple power-law (e.g.. iPerlman et al. 20051 ). in- 
dependent of flux level and - as demonstrated in this work - extending into the hard X-ray 
range. If interpreted as signatures of the particle energy gain process, one may conclude 
that acceleration takes place at all flux levels observed so far. In all cases there were no 
indications for a significant amount of X-ray absorbing material in the line-of-sight beyond 
the Galactic absorption column density. 

The average spectrum during our Suzaku observations was qui te similar to the on e 



from the 20% higher flux XMM observations in 2001 (r ~ 2.1 - 2.4: IPerlman et al. 20051 ) 



where also no sign ificant variability was noticed. The 1997/98 BeppoSAX observations 



(jWolter et al. 20001 ) had a similar value of the break energy at a factor ~ 1.3 — 1.9 higher 
flux level. The flux reported from the 2004 observations with XMM was comparable to the 
one measured by BeppoSAX in 1997, however, the spectral break shifted to lower energies 
by a factor ~ 1.2. Because of the lacking soft X-ray sensitivity, the PCA coverage of 1ES 
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1101-232 in 2005 revealed a break energy at much higher energies (~ 8 keV). If the 1- 
8 keV X-ray spectral index is compared to the observed flux level at each observation, 
no statistically convincing indication for a flux-index correlation is found (Spearman rank 
correlation coefficient R s = —0.34 with chance probability P c = 0.33, Kendall's tr = —0.30 
with P c = 0.23). Comparing, e.g., the XMM measurements carried out in 2001 and 2004, no 
significant spectral changes within the uncertainties occured despite a factor ~ 1.6 higher 
flux in 2004. On the other side, the BeppoSAX observations revealed a spectral softening 
from the higher flux state in 1997 to the lower flux state in 1998. When compared to the 
high flux RXTE data from 2005 a noticable softening occured despite the overall high flux 
level. If following the evolution of the spectrum versus flux, no clear spectral hysteresis 
could be found. In summary, it seems that spectral changes observed from annual visits of 
this source do not follow particular patterns or relations with respect to flux changes. The 
overall spectral changes lie within a rather narrow range between ~ 2 and 2.5 for the 1-8 keV 
photon spectral index. 

The measured so far comparably long X-ray variability time scale in 1ES 1101-232 may 
imply that this source was observed only in the quiescence over the last ~ 15 years. If 
one assumes the measured activity states from one year to the other to be not causally 
related, the annual measurements may then be treated as an independent ensemble of mea- 
surements from an HBL. In particular, it seems reasonable to probe the behaviour in the 
peak luminosity-peak energy diagram within the available HBL range. For a robust deter- 
mination of the peak energy, low energy data, preferably in the IR/optical/UV range, are 
necessary. Giommi et al. (2005) used optical data from the literature in combination with 
the BeppoSAX measurements, while both XMM pointings had the advantage of simultane- 
ous optical coverage with the Optical Monitor (OM) onboard the satellite, and SWIFT the 
onboard UVOT instrument. We did not find any convincing anti-correlation between peak 
energy and peak flux within the available 6 ensemble members (R s = —0.6 with P c = 0.21, 
r K = -0.33 with P c = 0.35). 

1ES 1553+113 is known as a generally X -ray soft, bright HBL with only modest X- 



ray variability (factor < 3.5). SWIFT in 2005 flTramacere et al. 20071 : iMassaro et al. 20081 ) 
achieved a detection of this source up to only 8 keV, despite its high flux level in October 
2005 that was reached w ithin half a year from a 3.5 times lower flux level. BeppoSAX in 



1998 (IDonato et al. 20051 ) saw it only with the LECS/MECS at a relatively low flux level, no 
detection with the PDS was possible. With our Suzaku observation, 1ES 1553+113's spec- 
trum could finally be measured up to 30 keV while descending from the 2005 high flux level. 
Like for 1ES 1101-232, past X-ray observations encompass mostl y snapshot observation s of 



typically 3-10 ksec each. RXTE's observation campaign in 2003 (lOsterman et al. 20061 ) in- 



volved visits to the source about 3 times per day for ~ 3 ksec, for a total period of 21 days. 
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During this time a very smooth rise in flux up to 3 times the lowest flux level was observed. 
No hints of subhour variations were recorded in any of the past observations. With our 
~ 50 ksec continuous Suzaku measurements we tested for intraday variability, but obtained 
a null result also at this intermediate flu x level which was quite comparable to the 2001 flux 



recorded by XMM flPerlman et al. 20051 ). 



No significant differences were found between the spectra reported from the 2001 and 
2006 observations: the average spectrum during our Suzaku observations was, within the un- 
certainties, in good a greement to the one from the XMM observations in 2001 (T ~ 2.2 — 2.4: 



Perlman et al. 20051 ). The surprising spe ctral stability of this s ource holds even for the low 



flux state recorded with RXTE in 2003 (jOsterman et al. 20061 ). despite the factor ~ 5 dif- 



ference in flux. Even within the tripling of the flux during the PCA observations spectral 
changes were not significantly beyond their uncertainties. In particular, any systematic spec- 
tral trend with changing flux is not apparent, neither on a year-by-year time scale nor within 
the rising part of the RXTE monitoring data within their admittedly large uncertainties: 
as is the case for 1ES 1101-232, the photon index above 3 keV typically lingered around 2 
and 2.5 even during flux increases by a factor ~ 2 or more. Those modest spectral changes 
together with the smoothness of the long-term light curve may suggest that explosive events 
are among the least likely scenarios to account for the observed flux variations here. 



The 1998 BeppoSAX observations (IDonato et al. 20051 ) had a lower break energy 



1 keV) at a factor ~ 2.5 lower flux level as compared to the 2006 Suzaku observations. 
SWIFT detected 1ES 1553+113 in both, a rather low flux state (April 2005) and relatively 
high flux level (October 200 5), but with a surpris ingly stable peak energy of ~ 0.4 — 0.5 keV 



in the XRT spectrum (e.g. iMassaro et al. 20081 ). ~ 3 times lower than the Suzaku data 



from 2006 imply. The flux reached during the April 2005 observations with SWIFT was 
comparable to the one measured by BeppoSAX in 1998, however the peak energy from the 
XRT data shifted to lower e nergies by a factor ^2 .5. No break energies were reported from 



the 2003 RXTE/PCA data flOsterman et al. 20061 ) . 



Following the argumentation above we probe here the dependence between peak lu- 
minosity and peak energy also for 1ES 1553+113, and utilize preferably data from the 
IR/optical/UV range to complement the X-ray band measurements. Among the past obser- 
vations only the XMM and SWIFT pointings were quasi- simultaneous measurements with 
optical/UV available, namely OM, UVOT and ROTSE data, respectively. We use here 
our Suzaku observations with the data from the quasi-simultaneous optical coverage by the 
KVA optical telescope for a combined X-ray - optical spectral fit. A log-parabolic shape 

was used to determine the peak energy to ~ 0.01 keV with parame- 
ters P = 0.075 ± 0.001 and T = 2.298 ± 0.003. Note, however, that this representation does 
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only give a Xrcd ~ 5.7 due to the large weight given to the only optical point (that has a small 
uncertainty) with respect to the many X-ray data points. More lower energy points would 
clarify this. We did not use the peak values given for the 2003 RXTE campaign for this 
undertaking to assure for a sample with equally determined peak energies: the cm radio data 
from Osterman et al. (2006) for determining the peak energy would give too much weight 
towards the longer wavelengths when compared to a procedure where only optical/UV in 
combination with the X-ray data are used. The peak energies from the UVOT measure- 
ments range from ~ 0.05 to 0.02 keV from the high to the low flux state, respectively. The 
latter is in very good agreement with the values derived from the XMM 2001 observations. 
Also for 1ES 1553+113 any convincing systematic correlation between peak energy and peak 
flux within the very limited sample of 4 ensemble members is lacking evidence (R s = 0.39 
with P c = 0.61, tk = 0.55 with P c = 0.26). This may indicate that particle acceleration 
up to the maximum energy is not preferentially limited by Compton losses if beaming does 
not change significantly. Indeed, because of the lack of sufficiently dense radiation fields in 
HBL-type AGN, synchrotron losses may dominate the radiative particle loss channel. The 
absence of a link between the peak luminosity and peak energy can then be understood if 
either an increase of the jet's synchrotron radiation field is not connected to an increase of 
the magnetic field strength alone on year time scales, or expansion losses and/or escape from 
the emission region sets a boundary to the maximum electron energy at the so far observed 
activity states in 1ES 1101-232 and 1ES 1553+113. 



5.2. Spectral energy distributions 

The overall SED of 1ES 1101-232, including the Suzaku spectrum obtained by us and 
the s imultaneous H.E.S.S. da ta is shown in Fig. O If corrected by the minimal EBL (P0.45 



from lAharonian et al. 2006al ). comparable power is emitted at 7-rays than at the synchrotron 
hump with a 7-ray peak beyond TeVs at the time of the quasi-simultaneous Suzaku - TeV 
observations. In the framework of a simple one-zone SSC model a (assumed uniform) field 
strength of sub-equipartition strength B 2 oc {yF v jc)/{yF^ syn ) can then be deduced. Al- 
ternatively, further target photon fields in the vicinity of the 7-ray emission zone can lead 
to an enhanced 7-ray photon output. In hadronic blazar emission models TeV-emission 
is the result of either reprocessed cascade emission initiated by pairs and/or photons pro- 
duced in photomeson int e ractions, or by pro ton and/or 7r ± //^-synchrotron radiation (e.g. 



Miicke fc Protheroe 200 ll ; iMiicke et al. 20031 ). In environments of weak ambient photon 



fields proton-photon interactions constitute a correspondingly sub-dominant contribution to 
the 7-ray component, which then leaves proton synchrotron (or curvature) radiation as the 
main 7-ray producer. 
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Fig. 9. — Broadband SED of 1ES 110 1-232. The observed H. E.S.S. data were de-absorbed 



using the minimal EBL (P0.45 from lAharonian et al. 2006al ). Present Suzaku data and 
the simultaneous de-absorbed July 2006 H.E.S.S. data are indicated in r ed. Blue points 
(XMM, RXTE, HESS) indicate data from 2004/05 jAharonian et al. 2007ah . while historical 
data collected from the literature are in black open symbo ls, and light grey (BeppoSAX: 
Wolter et al. 20001 ) . The SWIFT data (jMassaro et al. 20081 . not shown here) are compatible 
with the XMM flux level. 
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TeV observations of 1ES 1101-232 were used previously (jAharonian et al. 2006al ) to 



set the most stringent l imits on the EBL at ~ 0.7 — 4/xm using spectral considerations 



(jAharonian et al. 2006al ). These limits could not be improved owing to the limited quality 



from the low flux TeV-data from 2006. 

The overall SED of 1ES 1553+113 with the contemporaneous Suzaku, VHE (H.E.S.S., 
MAGIC) and optical (KVA) data is shown in Fig. 110} complemented with historical ones from 
the literature. Using a minimal EBL (P0.45) for de- absorption, approximately equal or more 
power is emitted at 7-rays than at the synchrotron hump with a 7-ray peak at sub-TeVs. In 
the framework of a simple one-zone SSC model, this ratio indicates sub-equipartition field 
strengths in the emission region. Furthermore, if the same electrons emit TeV- and X-ray 
photons co-spatially, then the intrinsic spectral index at VHEs must not be harder than 
at X-ray energies. For an EBL density at the galaxy counts flux level this is the case for 
z > 0.4. Alternatively, external photon fields may contribute to the target photon fields for 
inverse Compton scattering, or hadronic emission models may be at work here. 



-27- 



-10 - 



-11 



12 - 



-13 



14 - 



"i 1 1 r 



-| 1 1 r 



-i 1 1 r 




9i 

& 

□ 



1ES 1553+113 



O 



J L 



J L 



- 47 



- 46 



N 

45 „ 



go 

44 <d 



43 §? 



- 42 



10 



15 20 
Log i/ [Hz] 



25 



Fig. 10. — Broadband non-simultaneous SED of 1ES 1553+113. The quasi-simultaneous 
data points (KVA, Suzaku, H.E.S.S., MAGIC) data are indicated a s (red) squares. The 
green/grey line respresents the highest flux level observed by SWIFT ( jMassaro et al. 20081 ) 
of this source. The 7-ray data have been de-absorbed for a source redshift of z = 0.3 using 
the minimal EBL (P0.45 from lAharonian et al. 2006al ). 
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6. Summary 

Continuous ~ 50 ksec Suzaku pointings towards each of the two distant TeV-blazars, 
1ES 1101-232 and 1ES 1553+113, were carried out in May and July 2006, respectively, with 
quasi-simultaneous coverage at VHEs (H.E.S.S., MAGIC). In this work we presented the 
Suzaku data set, and set them into context to the VHE observations and to their behaviour 
in the past. Our findings can be summarized as follows: 

• Suzaku observations in 2006 of 1ES 1101-232 and 1ES 1553+113 found both TeV- 
blazars in a non-active state. We detected both blazars from soft to hard X-rays up 
to ~ 30 keV, with 1ES 1553+113 measured to so far unprecedented energy in the 
synchrotron emission component. 

• The observed X-ray flux of 1ES 1101-232 is the lowest ever measured from this source. 

• The combined XIS and PIN spectra of both sources showed surprisingly similar spectral 
shapes. This may indicate a common physical mechanism as being responsible for such 
spectral signature. Both show indications of curvature, and can be equally well fitted 
by a broken power law or log-parabolic shape. This curvature is found to extend 
well into the hard X-ray band (< 30 keV). Note, however, that the PIN data of 1ES 
1553+113 indicates a steepening beyond the log-parabolic or broken power law model. 

• Simultaneous H.E.S.S. and MAGIC observations during the Suzaku epoch revealed no 
variability on any time scale probed by our data. H.E.S.S. obtained a weak detection 
(~ 3a) of 1ES 1101-232 at the low flux level observed also during previous years. The 
combined 1ES 1553+113 VHE spectrum (90 — 500 GeV) did not show any significant 
changes with respect to historical VHE observations. Consequently, this data set did 
not allow to put further constraints on the EBL flux density. 

• 1ES 1553+113 is unlikely to reside at a redshift z > 0.7 even if the EBL is at the 
minimum flux density (inferred from integral galaxy counts). 

1ES 1553+113 and 1ES 1101-232, 7-ray blazars with surprisingly mild variability at 
virtually all energies, but nonetheless produce photons up to the highest TeV energies. This 
can be contrasted to the prominent TeV-blazars like Mkn 501, Mkn 421 or PKS 2155-304 
where rapid variations on subhour time scales at X- and TeV-energies have been detected at 
preferentially high activity states. Are there consequently various modes of high energy pro- 
duction in blazar jets realized in nature - rapidly explosive and quietly smooth ones - , or are 
we being mislead by observations at the instrument sensitivity limit and/or sparse sampling 
modes? Future high-sensitivity, dense sampling observations may answer this question. 
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